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This article discusses the way the standard description of capillary filling dynamics has to be modified to 
account for liquid/solid slip in nanometric pores. It focuses in particular on the case of a large slip length 
compared to the pore size. It is shown that the liquid viscosity does not play a role, and that the flow is 
only controlled by the friction coefficient of the liquid at the wall. Moreover in the Washburn regime, the 
filling velocity does not depend on the tube radius. Finally, molecular dynamics simulations suggest that this 
standard description fails to describe the early stage of capillary filling of carbon nanotubes by water, since 
viscous dissipation at the tube entrance must be taken into account. 



I. INTRODUCTION 

Dynamics of capillary rise is a standard case of cou- 
pling between capillarity and hydrodynamics^'^. The fill- 
ing behavior has been investigated since the beginning of 
the 20th century-^"^. In particular, the famous (Lucas- 
)Washburn law (neglecting the liquid inertia) states that 
the filling velocity decreases as the inverse square root of 
time. But this old problem needs to be revisited as re- 
gards nanometric pores^'^. Even if continuum hydrody- 
namics remains valid for simple liquids (including water) 
down to channel sizes of a typical Inm^*^, surfaces will 
play an increasing role. In particular, deviations from the 
classical hypothesis of a no-slip boundary condition (BC) 
at the liquid/solid interface have been predicted theoret- 
ically and observed experimentally^^. The simplest way 
to account for liquid/solid slip is the so-called 'partial 
slip' BC, which links the slip velocity iigUp with the shear 
rate at the solid surface dnv: tigiip = bdnV, in which the 
slip length b is the depth inside the solid where the linear 
extrapolation of the velocity profile vanishes (Fig. l.a). 
Let us emphasize here that the slip length, though it has 
a simple geometric meaning, is not a fundamental prop- 
erty of the liquid/solid interface. Indeed, the partial slip 
BC physically stems from the continuity of tangential 
stress at the wall: the viscous shear stress exerted by the 
liquid on the wall rj dnV {rj being the shear viscosity of the 
liquid) is equal to the friction force suffered by the liquid 
from the wall, which can be written as F/A = — Ai'siip 
where A is the liquid/solid friction coefficient and A the 
contact area. The slip length is accordingly a combina- 
tion of the bulk liquid viscosity and the interfacial fric- 
tion coefficient: b = rj/\. For simple liquids on smooth 
surfaces, slip lengths up to a few tens of nanometers 
have been experimentally measured^ ^. Liquid-solid slip is 
therefore expected to significantly affect flows in channels 
of nanometric size, even when continuum hydrodynamics 
remains valid. 
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In recent years, many works have investigated capillary 
filling at the nanoscale, emphasizing the roles of dynamic 
contact angle, liquid inertia, and liquid-solid slip^^^^^. 
In this context, recent experiments^^"^° and numerical 
simulations'^"'^''^^ have reported slip lengths of water (and 
other liquids) in carbon nanotubes (CNT) much larger 
than the tube radius. This article seeks to investigate the 
way this quite special condition modifies the dynamics 
of capillary filling. Firstly, the full equations describing 
capillary filling with a partial slip boundary condition 
will be derived and solved. The article will then focus on 
the case of CNTs where giant liquid/solid slip has been 
reported. Finally this standard analytical model will be 
tested against molecular dynamics simulations. 



II. ANALYTICAL MODEL 

Before turning to the specific case of water filling a 
CNT, the equation describing capillary filling with a par- 
tial slip boundary condition, taking into account the liq- 
uid inertia, will be derived and solved. This model is 
formally equivalent to the one introduced by Supple and 
Quirke"'^^, with the difference that the equation will here 
be derived in the framework of continuum hydrodynam- 
ics, as we have discussed that liquid/solid slip could be 
significant for channel sizes where continuum hydrody- 
namics remained valid. 

Let's consider a cylindrical pore of radius R, in contact 
with a reservoir of liquid (density p, viscosity rj), as illus- 
trated in Fig. l.a. The liquid surface tension is denoted 
7, and its contact angle on the pore wall 9 (it is assumed 
that the liquid is wetting the pore, i.e. < 7r/2). The 
distance of the meniscus from the pore entrance is de- 
noted L. Therefore the velocity of the meniscus V , equal 
to the average liquid velocity inside the pore, is related 
to L through: V = ^- The filling dynamics results from 
the competition between a capillary force at the contact 
line Fc, driving the motion of the liquid inside the pore, 
and a viscous friction force at the wall F^, slowing the 
liquid motion (gravity can usually be neglected at small 
scales). One can write the equation of motion for the 
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and; 



FIG. 1. (a) Schematics of the system considered, (b) Snap- 
shot of a typical system used for the Molecular Dynamics sim- 
ulations, created by using the VMD molecular visualization 
program . 



liquid inside the pore: 

djMV) 
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where M — pirR^L is the mass of liquid inside the pore. 
We will denote A7 — 730 ~ 7sl = 7 cos 9 the capillary 
force per unit length at the contact line, where 730 and 
73L are the solid-gas and solid-liquid surface tensions. 
The total capillary force is then simply given by: 



Fc = 2nRA'y. 



(2) 



The friction force is given by the Poiseuille law^'^, mod- 
ified to take into account liquid/solid slip'^'i^'^i.^s. 
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The equation of motion (1) can finally be re- written: 
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This non-linear equation admits the following solution-'^^: 
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One can easily check that in the long time limit, where 
the liquid inertia can be neglected compared to the vis- 
cous friction inside the tube, Eq. (5) simplifies into the 
well-known Washburn law'''^*'^^'^^: 
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On the contrary, in the short time limit viscous friction 
can be neglected compared to inertia. Eq. (5) simplifies 
into: 
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Of course neither the viscosity 77 nor the slip length 6 ap- 
pear here since the friction term has been neglected. One 
could finally note that even this inertial regime does not 
appear instantaneously, essentially because the meniscus 
shape changes when the liquid enters into the pore. But 
it can be shown that the meniscus adopts a stationary 
shape for L ~ R/2 (from^). Therefore, the pre-inertial 
regime disappears extremely fast in nanometric pores. 

Water inside CNT- Recent experiments^^"'^^ and 
numerical simulations'^^''^^ have reported surprisingly 
large slip lengths of water (and other liquids) in car- 
bon nanotubes. The first experiments^^ indicated slip 
lengths up to tens of micrometers, even if more recent 
expcrimcntaP^ and numcricaP^ works point to a smaller 
effect, with 6 on the order of a few hundredth of nanome- 
ters. When the slip length is much larger than the pore 
size: b ^ R, the flow inside the pore will be almost 
plug-like, namely with a constant velocity profile. In this 
regime, there is almost no shear inside the liquid, and 
the viscosity is not expected to play a role. In fact, the 
viscous term F^ will only depend on the intcrfacial fric- 
tion coefficient A. This can be seen if one take the limit 
b/R > 1 in Eq. (3): 
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But this expression can be found directly under the 
assumption of a perfect plug flow, where the velocity 
of the liquid is everywhere equal to its average value 
V. The friction force over the liquid/solid contact area 
A = 2'kRL is then simply: 
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Generally the ratio ?7/(l+46/i?) can be replaced by Ai?/4 
in the equation of motion (1), its solutions Eqs. (8) and 
(5), and the expression of the transition time (6) and 
transition length (7). In particular the Washburn law (in 
the viscous regime) becomes: 
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In the plug-flow case, viscosity does not play a role, and 
the flow is only controlled by the friction coefficient. Fur- 
thermore, the filling dynamics does not depend on the 
tube radius anymore. This important difference could be 
used as a clear experimental signature of a plug flow in 
the viscous regime. 

To compute orders of magnitude, one also needs to 
quantify the driving force Fc- This is far from triv- 
ial: Whereas measurements of water contact angle on 
highly oriented pyrolytic graphite converge to a value 
of 86 °'^'*'"^^, corresponding to 6* = 96° for water on a 
single graphene sheet (hence a non-wetting situation, 
where water should not invade CNT), there are exper- 
imental evidences that water fills at least the very small 
CNTs^^''^^"'^^. At these scales however, it is hard to apply 
macroscopic concepts of capillarity like surface tension 
and contact angle"'". But if water indeed fills CNT, Fc 
will be positive and one can estimate that it will be at 
most on the order of 27ri?7. 




B=0.514nm 

— fl= 0.853 nm 

— fl= 1.19 nm 

— fl= 1.87nm 
-- linear fits 



0.1 0.2 0.3 

tins) 



-I — I I I I I 1 1 — r 



© ® 



C^l I I I I I I I I I L 



R(nm) 



FIG. 2. (a) Evolution of tlie filling length L with time t, for 
different tube radii. Dashed lines represent two examples of 
linear fit used to compute the filling velocity V. (b). Filling 
velocity 1/ as a function of tube radius R: Simulation data 
(circles) and prediction of the standard model (line), using 
Eq. (9). 



III. MOLECULAR DYNAMICS SIMULATIONS AND 
DISCUSSION 

We considered a water reservoir in contact with an 
initially empty CNT (Fig. l.b). The CNT length was 
10 nm, with radii ranging between 0.514 and 1.87 nm. 
The reservoir was bordered with two graphene sheets, 
at a distance large enough to ensure that the liquid wa- 
ter was always in equilibrium with its vapor. Periodic 
boundary conditions in the two directions perpendicular 
to the tube axis ensured that the water surface remained 
planar. Consequently, the pressure in the water reservoir 
always remained constant, at a value extremely close to 
zero. 

The Amber96 force field^^ was used, with TIP3P wa- 
ter and water-carbon interaction modeled by a Lennard- 
Jones potential between oxygen and carbon atoms, with 
parameters £oc = 0.114kcal/mol and aoc = 3.28 A. The 
tabulated density, surface tension and viscosity of TIP3P 
water at 300 K and liquid/vapor coexistence are respec- 
tively: p = 0.980g/cm3, 7 = 0.0523N/m (from^^)^ ^nd 
77 = 0.321 mPas (from^"^). The contact angle of water on 
a single graphene sheet was measured independently to 
he = 57° for this model. The simulations were car- 
ried out using LAMMPS^"'. Long-range Coulomb forces 
were computed using the particle-particle particle-mesh 
(PPPM) method; A timestep of 2 fs was used, unless 
specified. The positions of the carbon atoms were fixed 
(simulations with fiexiblc and fixed walls were shown to 
give similar results for the statics and friction of confined 
liquids'^^'"^^'"'^). Water molecules were kept at a constant 
temperature of 300 K using a Nose- Hoover thermostat, 
applied only to the degrees of freedom perpendicular to 
the tube axis, with a damping time of 200 fs. Alterna- 
tively a Dissipativc Particle Dynamics thermostat was 
used (see details in the following). 

Water molecules are initially disposed on a simple cu- 



bic lattice with equilibrium density, so firstly the system 
is equilibrated during typically 120 ps, with a plug at 
the tube entrance to prevent water from entering. Then 
the plug is removed and the evolution of the water mass 
inside the tube is recorded as a function of time. To com- 
pute the filling length, the mass is divided by the linear 
density of water inside the tube, measured once the tube 
is completely filled. Results for different tube radii are 
presented in Fig. 2. a. It appears that L scales linearly 
with time, as expected from the model in the short-time 
limit. In Fig. 2.b, the filling velocities (slope of the L{t) 
curves) for different radii are compared to the theoreti- 
cal prediction of Eq. (9). The theoretical prediction fails 
to describe the numerical data both quantitatively (sim- 
ulated velocities are much smaller than expected) and 
qualitatively: While the model predicts that V should 
decrease as i?~^/^, the simulation shows that V does not 
depend significantly on the tube radius. 

While the quantitative discrepancy could be explained 
by a reduced surface tension (due to curvature and in- 
teraction with the surface), the invariance of V with the 
radius R seems to point toward another effect, for in- 
stance a drag force missing in the model. Even if vis- 
cosity is not expected to play a role inside the tube due 
to the very large slip length, its influence was neverthe- 
less tested, using a Dissipative Particle Dynamics (DPD) 
thermostat*^. This amounts to adding pairwise interac- 
tions between atoms, with a dissipative force depending 
on the relative velocity between each pair and a random 
force with a Gaussian statistics. This method has the 
advantage of preserving hydrodynamics. Furthermore, 
the amplitude of the dissipative term can be tuned to 
modify the liquid viscosity without changing its static 
properties. The viscosity dependency on the dissipative 
term was carefully calibrated in independent simulations 
of shear flow of water between two silica surfaces (Fig. 
3. a). For the largest dissipative terms, the timestep had 
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FIG. 3. (a) Calibration of viscosity 77 versus DPD coefficient 
7: MD results (circles) and fitting curve (line), see text for 
detail, (b) Filling velocity V as a, function of liquid viscosity 
77: MD results (circles) and fit with an inverse function (line). 



to be reduced down to 0.5 fs. 

The filling velocity was then measured for a given 
CNT, with radius 1.87 nm, varying the viscosity of the 
DPD-thermostated-water. On Fig. 3.b, one can observe 
a clear influence of viscosity on the filling velocity. Due 
to the very large slip length of water inside CNT, friction 
at the wall should be completely negligible for the filling 
lengths simulated. Moreover, in such a plug-flow regime, 
the liquid moves as a block everywhere inside the tube. 
Therefore, there should be no dissipation near the con- 
tact line, and consequently no dynamical change of the 
contact angle. Yet, even in the absence of friction inside 
the tube, viscous dissipation can occur at the tube en- 
trance due to the contraction of the streamlines^^'^^"^^. 
The power dissipated by viscosity can be estimated as: 
Pv = 2 /v '^"^(^^^J "'" '^j'^if- The only length scale re- 
lated to the entrance of the liquid inside the tube is the 
tube radius, hence: Pv « riR?{V/RY oc rjRV^. Ne- 
glecting the kinetic energy of the liquid (inertial term): 
Pk oc pB?V^ , the viscous dissipation must be equal to 
the power of the driving capillary force: Pc = 2tiR/S.'^V . 
Finally P^ = P^ ^ V (x A7/77. Therefore, if the fill- 
ing is controlled by viscous dissipation at the entrance, 
the velocity should not depend on the tube radius (as 
observed in Fig. 2.b), and scale with the inverse of vis- 
cosity. One can check in Fig. 3 that the numerical results 
for V versus 77 are correctly fitted with an inverse func- 
tion. To conclude, MD simulations show that, for water 
filling CNTs, inertia is negligible with regard to viscous 
dissipation at the entrance. The standard description of 
capillary filling presented earlier must thus be modified. 

In this paragraph, a simple analytical description of 
the full capillary filling dynamics, taking into account 
entrance effects (but neglecting liquid inertia), will be 
presented. It will be restricted to the plug-flow limit (6 » 
P), which is relevant to the case of water inside CNTs, 
but could easily be generalized. The pressure drop along 
the liquid inside the tube is given by the Laplace pressure 
difference across the curved liquid/vapor meniscus, minus 



the pressure drop at the tube entrance"''^: Ap = 2A7/P— 
StttiV/R. If we neglect liquid inertia, the corresponding 
velocity in the plug-flow regime is: V = RAp/(2XL). 
The resulting equation: 



2XL—- + imi-r- = 2A7 
at at 

can be integrated into: 

\L^ + Sttt^P = 2A7 t. 



(13) 
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This simple quadratic equation admits the following so- 
lution: 



where: 



and: 
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Lc = Stttj/X = Sirb. 



(17) 

For t -^ Tc, the full solution (15) simplifies to give a filling 
at constant velocity V = 2A^/{3TTr]), which is compatible 
with the simple dimensional analysis presented above. 
For i ^ Tc it reduces to Eq. (12). Interestingly, even in 
the plug-flow limit, the transition length between the two 
regimes provides a direct measure of the slip length (Eq. 
17). Furthermore, for slip lengths of a few hundredth 
of nanometers^®''^^, Lc would lie within the micrometer 
range. 



IV. CONCLUSION 

In this article, the consequences of the giant slip of wa- 
ter inside carbon nanotubes on the dynamics of capillary 
flUing have been investigated. It has been shown that 
in case the slip length is much larger than the tube ra- 
dius (plug-flow limit), viscosity does not play a role, and 
the relevant physical parameter controlling the flow is 
the liquid/solid friction coefficient. Moreover, scaling of 
the various quantities as a function of the tube radius are 
modified. It was then shown, using MD simulations, that 
in the short-time limit the filling velocity is not limited 
by the liquid inertia, but rather by viscous dissipation 
at the tube entrance. A model of the full capillary fill- 
ing dynamics, taking into account entrance effects (but 
neglecting liquid inertia) was finally presented. 

Generally, friction of water inside CNT is so small that 
for any real flow experiments, one can expect that the 
hydrodynamic resistance at the tube input and output 
will play a significant role^^. It could even contribute 
much more strongly to the total pressure drop than the 
friction inside the tube. In any case, measurement of 
flow rate versus pressure drop for CNTs of a given length 
can only give a higher limit for the friction coefficient 
inside the tube, so that measurements with different tube 
lengths could prove necessary. 
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